AD-A285  586 

"  •  .  _  _ _ .....  BHI  IkSI 


NUWONPT  Technical  Report  10,719 
22  August  1994 


A  Model  of  an  Infinite  Length, 
Axisymmetric,  Isotropic, 
Forced  Thick  Sheli 

Andrew  J.  Huli 
Submarine  Sonar  Department 


OGT  1  Gl934j^ 


Naval  Undersea  Warfare  Center  Division 

Newport,  Rhode  Island 


Approved  for  public  release;  distribution  is  unlimited. 


PREFACE 


This  investigation  described  in  this  report  was  sponsored  by  the  Office  of  Naval 
Research. 

The  technical  reviewer  for  this  report  was  S.  A.  Austin  (Code  2141). 

The  author  wishes  to  thank  Karen  Holt  (Code  0251),  NUWC  Detachment  New  London, 
for  her  help  with  the  technical  editing. 


REVIEWED  AND  APPROVED:  22  August  1994 


J.  Martin 

Acting  Head,  Submarine  Sonar  Department 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMBNo.  0704-0188 


MtiirMad  lo  avarad*  1  hour  par  raaponta.  including  lha  lima  lor  ravieiMng  inetrudians.  saarchmg  anating  data  aourcas,  galhanng  and 
mainlalning  tha  data  noadad.  and  comptating  and  raviewing  lha  cotoclion  «  Mormation.  Sand  cormnants  lagarding  ths  burdan  estimala  or  any  olhar  aspaci  of  ihB  ooMaction  ol  intormation. 
including  tuggaatiana  lor  raducing  this  burdan,loWaahlnglanHaadquanaraSarvicaa,Oiraclarala  lor  intormatlonOparalions  and  Rapona.  1215  Jallarson  Oavia  Highway.  Suae  1204.  Arlington. 

and  to  tha  Ollioe  ol  Manaoemant  and  Budoat.  Paparwork  Raduciion  Protect  f0704-0l88i.  Washirioion.  DC  20503. _ 


AGENCY  USE  ONLY  (Leave  Blank) 


2.  REPORT  DATE 

3.  REPORT  TYPE  AND  DATES  COVERED 

22  August  1994 

Final 

4.  TITLE  AND  SUBTITLE 

A  Model  of  an  Infinite  Length,  Axisymmetric,  Isotropic, 
Forced  Thick  Shell 


6.  AUTHOR(S) 

Andrew  J.  Hull 


PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Undersea  Warfare  Center 
Detachment  New  London 
New  London,  Connecticut  06320 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

TR  10.719 


SPONSORING/MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

Office  of  Naval  Research 
800  North  Quincy  St. 

Arlington,  VA  22217 


10.  SPONSORING/MONITORING  AGENCY 
REPORT  NUMBER 


12a.  OISTRIBUTION/AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

The  model  described  in  this  report  provides  an  exact  solution  to  the  equations  of  motion  for  an  infinite  length, 
axisymmetric,  isotropic,  forced  thick  shell.  In  this  model,  the  governing  equations  of  motion  of  an  isotropic  elastic 
solid  in  cylindrical  coordinates  are  first  separated  into  a  dilatational  wave  equation  and  three  distortional  wave 
equations  that  can  be  solved  with  Bessel  functions.  These  general  solutions  are  then  inserted  into  stress-strain 
constitutive  equations  at  the  interior  and  exterior  of  the  shell.  The  displacements  are  next  determined  from  forcing 
functions  that  are  applied  to  the  shell.  In  this  new  model,  the  shell  thickness  effects  on  wave  propagation  are 
addressed.  The  case  where  the  mean  radius  of  the  shell  is  equal  to  its  thickness  has  been  analyzed  and  compared 
with  thin  shell  equations  of  motion.  It  is  shown  that,  in  general,  the  thin  shell  equations  of  motion  often  uncerpredict 
thick  shell  response. 


14.  SUBJECT  TERMS 

Axisymmetric  Shell,  Elastic  Solid,  Equations  of  Motion,  Forced  Thick  Shell,  Infinite  Length 
Shell,  Isotropic  Solid,  Stress-Strain  Equations,  Thick  Shell,  Thin  Shell,  Wave  Equation 


is.  NUMBER  OF  PAGES 

44 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CUSSIFICATION  OF  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  ABSTRACT 


UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


TR  10,719 


TABLE  OF  CONTENTS 

Page 


LIST  OF  ILLUSTRATIONS .  ii 

1  INTRODUCTION .  1 

2  SYSTEM  MODEL . 2 

3  MODEL  VALIDATION .  15 

4  ANALYSIS .  20 

5  CONCLUSIONS .  30 

6  REFERENCES .  30 

APPENDIX  -  DERIVATIVES  OF  THE  BESSEL  FUNCTIONS .  A-1 


1 


TR  10,719 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1  Cylindrical  Shell .  3 

2  Relationship  of  Bessel  Function  Regions  to  Wave  Speeds .  7 


3  Magnitude  of  Axial  Displacement  Divided  by  Applied  Axial  Stress: 

Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/=  50  Hz,  hla  =  1/100 .  18 

4  Magnitude  of  Radial  Displacement  Divided  by  Applied  Axial  Stress: 

Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/=  50  Hz,  hla  =  1/100 .  18 

5  Magnitude  of  Axial  Displacement  Divided  by  Applied  Radial  Pressure: 

Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/=  50  Hz,  hJa  =  1/100 .  19 

6  Magnitude  of  Radial  Displacement  Divided  by  Applied  Radial  Pressure: 

Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/=  50  Hz,  hla  =  1/100 .  19 

7  Magnitude  of  Axial  Displacement  Divided  by  Applied  Axial  Stress: 

Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/  =  50  Hz,  2hl{a+b)  =  1...  22 

8  Magnitude  of  Radial  Displacement  Divided  by  AppUed  Axial  Stress: 

Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/=  50  Hz,  2hl{a-¥b)  =  1...  22 

9  Magnitude  of  Axial  Displacement  Divided  by  Applied  Radial  Pressure: 

Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/=  50  Hz,  2hl{a+b)  =  1...  23 

1 0  Magnitude  of  Radial  Displacement  Divided  by  Applied  Radial  Pressure: 

Comparison  of  Thick  Shell  Model  to  TTiin  Shell  Model,/=  50  Hz,  2hl{a+b)  =  1...  23 

1 1  Magnitude  of  Axial  Displacement  Divided  by  Applied  Axial  Stress: 

Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model, /=  500  Hz,  2hlia+b)  =  1..  24 

1 2  Magnitude  of  Radial  Displacement  Divided  by  Applied  Axial  Stress: 

Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/=  500  Hz,  2hlia+b)  =  L.  24 

13  Magnitude  of  Axial  Displacement  Divided  by  Applied  Radial  Pressure: 

Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/=  500  Hz,  2}il(a+b)  =  1..  25 

14  Magnitude  of  Radial  Displacement  Divided  by  Applied  Radial  Pressure: 

Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/=  500  Hz,  2hl{a+b)  =  1..  25 

1 5  Magnitude  of  Axial  Displacement  Divided  by  Applied  Axial  Stress: 

Comparison  of  Transfer  Function  at  r  =  o  to  r  =  b,/=  50  Hz,  2hlia+b)  =  1 . 26 


ii 


TR  10,719 


LIST  OF  ILLUSTRATIONS  (CONT'D) 

Figure  Page 

1 6  Magnitude  of  Radial  Displacement  Divided  by  Applied  Axial  Stress: 

Comparison  of  Transfer  Function  at  r  =  a  to  r  =  f>,/=  50  Hz,  2h/ia+b)  =  1 . 26 

17  Magnitude  of  Axial  Displacement  Divided  by  Applied  Radial  Pressure: 

Comparison  of  Transfer  Function  at  r  =  a  to  r  =  fc,/=  50  Hz,  2hlia+b)  =  1 . 27 

1 8  Magnitude  of  Radial  Displacement  Divided  by  Applied  Radial  Pressure: 

Comparison  of  Transfer  Function  at  r  =  o  to  r  =  b,f=  50  Hz,  2hl{a+b)  =  1 . 27 

19  Magnitude  of  Axial  Displacement  Divided  by  Applied  Axial  Stress: 

Comparison  of  Transfer  Function  at  r  =  a  to  r  =  b,f=  500  Hz,  2h/(a+b)  =  1 .  28 

20  Magnitude  of  Radial  Displacement  Divided  by  Applied  Axial  Suess: 

Comparison  of  Transfer  Function  at  r  =  a  to  r  =  b,f=  500  Hz,  2hJia+b)  =  1 .  28 

2 1  Magnitude  of  Axial  Displacement  Divided  by  Applied  Radial  Pressure: 

Comparison  of  Transfer  Function  at  r  =  a  to  r  =  b,f=  500  Hz,  2hJ(a+b)  =  1 .  29 

22  Magnitude  of  Radial  Displacement  Divided  by  Applied  Radial  Pressure: 

Comparison  of  Transfer  Function  at  r  =  a  to  r  =  b,f  =  500  Hz,  2h/{a+b)  -  1 .  29 


iii/iv 

Reverse  Blank 


TR  10,719 


A  MODEL  OF  AN  INFINITE  LENGTH,  AXISYMMETRIC, 
ISOTROPIC,  FORCED  THICK  SHELL 

1.  INTRODUCTION 

The  displacement  field  and  wave  propagation  of  elastic  shells  have  been  studied  by 
numerous  researchers  because  of  the  use  of  such  shells  in  many  mechanical  designs.  A 
common  solution  is  to  model  the  shell  as  unforced  and  infinitely  long  in  order  to  determine  its 
free  wave  propagation  characteristics  (Herrmann  and  Mirsky,  1956;  Mirsky,  1964a).  Another 
technique  approximates  the  shell  equations  of  motion  and  derives  a  forced  response  (Mirsky, 
1964b;  McNiven  et  al.,  1966).  Although  the  free  wave  propagation  approach  is  useful  for 
unforced  systems,  it  does  not  account  for  shell  loading  at  the  boundary.  However,  the  forced 
wave  method  is  based  on  approximations  to  the  system  that  are  not  always  accurate. 

In  this  report,  the  analysis  of  an  infinite  length,  axisymmetric,  isotropic,  forced  thick 
shell  is  presented.  The  derivation  begins  with  the  partial  differential  equations  of  motion  of  an 
isotropic  elastic  solid  in  cylindrical  coordinates.  The  shell  displacements  are  written  as  a  sum 
of  a  dilatational  component  and  an  equivoluminal  vector  component  This  approach  allows  the 
equations  to  be  separated  into  a  dilatational  wave  equation  and  three  distortional  (shear)  wave 
equations.  General  solutions  to  these  four  wave  equations  are  then  determined.  These  general 
solutions  are  inserted  into  the  stress-strain  relationships  at  the  inner  and  outer  diameter  of  the 
shell  and  equated  to  the  externally  applied  loading  function  on  the  shell.  Based  on  these 
relationships,  a  four-by-four  system  of  linear  equations  is  developed.  The  solution  to  the  linear 
equations  yields  the  displacements  of  the  shell. 

Corresponding  to  an  extremely  long  cylinder  with  external  forcing  at  a  definite 
wavenumber  and  frequency,  the  model  has  been  derived  to  address  shell  thickness  effects  on 
wave  propagation.  Two  specific  cases  are  investigated;  the  first  is  an  extremely  thin  shell  and 
the  second  is  a  shell  where  the  mean  radius  is  equal  to  the  thickness. 


1 


TR  10,719 


2.  SYSTEM  MODEL 

The  system  model  is  a  cylindrical,  linear,  isotropic,  elastic  medium  whose  motion  is 
governed  by  the  equation  (Timoshenko  and  Goodier,  1934) 

/iV2u  +  (A+//)VV«u  =  p^  ,  (1) 

dr 

where  p  is  the  density;  A  and  p  are  the  Lam^  constants;  t  is  time;  •  denotes  a  vector  dot 
product;  u  is  the  cylindrical  coordinate  displacement  vector  expressed  as 

Ur(r,9,z,t) 

u  =  - UQir,d,z,t)  ■  ,  (2) 

Uyir,6,z,t) 

with  subscript  r  denoting  the  radial  direction,  6  denoting  the  angular  direction,  and  z  denoting 
the  axial  direction;  V  is  the  gradient  vector  differential  operator  written  in  cylindrical 
coordinates  as  (Potter,  1978) 

with  if.  denoting  the  unit  vector  in  the  r  direction,  ig  denoting  the  unit  vector  in  the  9  direction, 
and  denoting  the  unit  vector  in  the  z  direction;  V  is  the  three-dimensional  Laplace  operator 
operating  on  vector  u  as 

with  V  operating  on  scalar  u  as 

and  the  term  V  •  u  is  called  the  divergence  and  is  equal  to 

_  dur  1  dug  du,  Uf  ft-. 

The  coordinate  system  of  the  shell  is  shown  in  figure  1. 


2 


3 


TR  10,719 

The  displacement  vector  u  is  written  as 

u  =  V0  +  VxH,  (7) 

where  0  is  a  dilatational  scalar  potential,  x  denotes  a  vector  cross-product,  and  H  is  an 
equivoluminal  vector  potential  expressed  as 


H  = 


Hrir,e,z,t) 

HQ(r,6,z,t) 

H^{r,d,z,t) 


(8) 


Expanding  equation  (7)  and  breaking  the  displacement  vector  into  its  individual  terms  yields 


_  d0  1  dH,  dHg 
^  dr  r  dd  dz  ' 


“0 


^\d^  dHr  dH, 
r  dd  dz  dr 


(9) 

(10) 


and 


and 


„  1  dHr 

^  dz  r  dr  r  de  ' 

Equation  (7)  is  inserted  into  equation  (1),  which  results  in 

dr 


c:?v2h  =  ^. 
"  dp- 


(11) 


(12) 


(13) 


The  constants  and  are  the  dilatational  and  ^hear  wave  speeds,  respectively,  and  are 
determined  by 


X  +  2u 
V  P 


(14) 


and 


-F 

'p 


Cj  —  „  I 


(15) 
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The  conditions  of  infinite  length,  axisymmetric  response  (n  =  0),  and  steady-state  response  are 
now  imposed,  allowing  the  scalar  and  vector  potential  to  be  written  as 

0  =  g{r)cos{nd)e‘^e‘^  =  ,  (16) 

Hr  =  hrir)sin{ne)e‘‘^e‘^  s  0  ,  (17) 

He  =  he{r)cosine)e^^e‘^  =  he(r)e^e^^  ,  (18) 


H,  =  /i2(r)sin(/i0)<>'^e'^  s  0  ,  (19) 

where  k  is  the  wavenumber  of  excitation,  (O  is  the  frequency  of  excitation,  and  /  is  the  square 
root  of  -1.  Note  that  for  axisymmetric  response,  the  equations  of  motion  are  dependent  only 
on  the  scalar  potential  0  and  angular  contribution  Hq  of  the  vector  potential.  Additionally, 
because  //^  =  0  and  H^  =  0,  and  Hq  aid  0  are  not  functions  of  0,  equation  (10)  becomes 

<-8  =  ^  =  0.  (20) 

where  ( )  denotes  any  function.  Inserting  equations  (16)-(19)  into  equations  (12)  and  (13) 
yields  the  following  four  wave  equations: 


dr^  r  dr  [c^  J’ 
dr^  r  dr  (cf  r^  J 
r  dr  [cf  r‘^  J 


d\y)  tdh,(r)  f:4(r)-0 

dr  f*  ’ 
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Equations  (21)  and  (23)  are  now  solved  with  Bessel  functions.  No  solution  is  found  to 
equations  (22)  and  (24)  because  they  do  not  contribute  to  the  axisymmetric  response.  The 
solution  to  equation  (21)  is 


CiJo(ar)  +  C2yo(ar), 


S(r)  = 


Ci/o(ar)  +  C2KQ(ar), 


(25) 


where  Jq  is  a  zero-order,  first-kind  standard  Bessel  function;  Tq  is  a  zero-order,  second-kind 
standard  Bessel  function;  Iq  is  a  zero-order,  first-kind  modified  Bessel  function;  Kq  is  a  zero- 
order,  second-kind  modified  Bessel  function;  C|  and  C2  are  constants  (determined  below);  and 


(26) 


The  solution  to  equation  (23)  is 


C3Ji(/3r)  +  C4ri(Pr), 


he(r)  = \ 


(27) 


C3/i(/3r)  +  C4/fi()3r),  |itl>  — . 

where  7|  is  a  first-order,  first-kind  standard  Bessel  function;  Tj  is  a  first-order,  second-kind 
standard  Bessel  function;  is  a  first-order,  first-kind  modified  Bessel  function;  is  a  first- 
order,  second-kind  modified  Bessel  function;  C3  and  C4  are  constants  (determined  below);  and 


(28) 


Based  on  the  solutions  shown  in  equations  (25)  and  (27),  three  distinct  wavenumber  regions 
are  dependent  on  the  dilatational  and  shear  wave  speeds.  In  these  low,  medium,  and  high 
wavenumber  regions,  the  solution  set  changes  between  standard  and  modified  Bessel  functions 
as  shown  in  figure  2. 
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1 

1 

hQ{r)=  \ 

1  !  1 

1  1 
•  Angular  Potential  hgCr)  • 

1 

1 

1  A0(r)  = 

C3h(Jir)+  i 

1  C3li(Pr)  + 

C4Ki(Jir)  1 

hQir)  =  C3Jiilir)  +  C4Yi(JJr) 

1  C^K^Uir) 

r  '  T 

__J _ J _ J _ J _ J _ ^ 

-CD  -CD  CD  CD 

—  —  k  =  0  —  — 

Wavenumber,  k 

Figure  2.  Relationship  of  Bessel  Function  Regions  to  Wave  Speeds 
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The  displacements  and  external  forces  are  now  equated  by  use  of  the  stress-strain 
constitutive  equations  on  the  forced  surfaces  of  the  shell.  The  normal  stress,  strain,  and  radial 
forces  acting  on  the  cylinder  are  related  by 

arr{a,d,z,t)-{^  +  2n)err{a,6,z,t)  +  XeQffia,  0.z,t)  +  {a,d,z,t)  =  p{a,6,z,t)  (29) 

and 

arr{b,d,z,t)-{^^2p )e„ ib,d,z,t)  +  Xeggib, d,z,t)  +  ib,d,z,t)  =  p{b, 0,z,t)  ,  (30) 

where  arf{a,b\0',zj)  is  the  normal  radial  stress,  £„{a,b\0,z,t)  is  the  normal  radial  strain, 
eQQ(a,h;0,z,t)  is  the  normal  circumferential  strain,  £^(a,b;0,z,t)  is  the  normal  longitudinal 
strain,  p{a,b,0,z,t)  is  the  external  pressure  on  the  shell  in  the  radial  direction,  and  a  and  b 
denote  the  inner  and  outer  radius,  respectively.  The  shear  stress,  strain,  and  longitudinal 
forces  are  related  using 

(a,  0,  z, t)  -  IjJLEfy  (fl,  0, ZJ)  =  0,z, t)  (31) 

and 

ar^ib,0,z,t)  =  2p£r,{b,0,z,t)  =  fib,0,z,t)  ,  (32) 

where  <7rz(a,b-,0,z,t)  is  the  shear  stress,  £,^{a,b;0,z,t)  is  the  shear  strain,  and  f{a,b;0,z,t) 
is  the  external  shear  stress  on  the  shell  in  the  longitudinal  direction. 

The  strains  are  related  to  the  displacements  in  an  axisymmetric  solid  by  (Timoshenko 
and  Goodier,  1934) 


dUr 

(33) 

(34) 

du. 

(35) 

and 
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'^~2[dr^  dzj‘ 


(36) 


The  relationship  between  the  displacements  (and  the  derivatives  of  the  displacements)  and  the 
potential  functions  g  and  hg  is  found  by  combining  equations  (9),  (11),  (16),  and  (18)  to 

produce 


ikZgitot 


_( d^gjr)  _  ..  dhQir)\  jtz  im 
dr  dr  f  ‘  • 


dz 

du^ 

dr 


^  1  dhgjr)  _  hgjr)  ^  d\y)\ 


dr  r  dr 


dr^ 


JkZei(ot 


(37) 

(38) 

(39) 

(40) 

(41) 


and 


\eikZeioM 


(42) 


Combining  equations  (29),  (30),  (33),  (34),  (35),  (37),  (39),  and  (40)  yields  the 
normal  stress  in  terms  of  the  potential  functions  g  and  hg  at  r  =  a  and  r  =  b.  These  equations 
are 


and 


dr  a  dr  dr 


(A  +  -  Kk^g(b)  -  , 

dr  b  dr  dr 


(43) 


(44) 


where  Pq  is  the  magnitude  of  the  normal  force  acting  on  the  interior  of  the  shell  and  is  the 
magnitude  of  the  normal  force  acting  on  the  exterior  of  the  shell.  Combining  equations  (31), 
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(32),  (36),  (41),  and  (42)  yields  the  shear  stress  in  terms  of  the  potential  functions  g  and  hg  at 


r  =  a  and  r-b.  The  results  are 


dr 


a  J  a  dr 


d^hgia) 


dhgib)  .  ..  d%[b)  _  „ 
dr 


where  Fg  is  the  magnitude  of  the  (externally  applied)  shear  stress  acting  on  the  interior  of  the 
shell  and  is  the  magnitude  of  the  (externally  applied)  shear  stress  acting  on  the  exterior  of 
the  shell.  Implicit  in  equations  (43)-(46)  is  the  assumption  that  the  external  loads  on  the  shell 
are  occurring  at  a  definite  frequency  and  wavenumber. 

Inserting  equations  (25),  (27),  and  (A'1)-(A-12)  from  the  appendix  into  equations  (43)- 
(46)  and  then  rewriting  as  a  four-by-four  system  of  linear  equations  results  in 

<*11  <*12  <*13  <*14  (Q  f^fl 

<*21  «22  <*23  <*24  Q  Pf, 

<*31  <*32  <*33  <*34  <-3 

<*41  a^2  <*43  <*443  IQJ  [fb. 

where  the  matrix  coefficients  are  determined  by  the  magnitude  of  the  wavenumber  k.  In 
the  low- wavenumber  region, 

!*!<-<-,  (48) 

and  the  matrix  coefficients  are 


a„=  (-A-2/i)  —  yo(«<*)  + 


(=r) 


/l(afl)-KA  +2/i)  Y  \J2{aa)  ,  (49) 


fll2=  (-A-2/t)  -y  -A/:^  yo(«<*)  + 


(^) 


ri(afl) 4-  (A  -H  2/i)  ^  y2(a<*)  >  (50) 


ai^  =  {-ilikp)Jo(J5a)  +  iinkl5)J2{pa)  , 
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ai4=H^ikP)YoiPa)  +  ii^kp)Y2{Pa)  , 


f  2  \ 

021=  (-A-2/i)  Y  7o(o*)+[^^j7l(afc)  +  (A  +  2/i)  ^Wab)  ,  (53) 


022=  (-A-2/i)  ^  -Xk^\YQ(ab)+ 


+  (A  +  2/i)  ^  l'2(afc)  .  (54) 


023  =  i-itikp)Jo(Pb)  +  mPVliPb)  ,  (55 

024  =  {-mkp)YQ{Pb)  +  mP)Y2(Pb)  ,  (56) 

031  =(-/2/i*a)yi(ao)  .  (57) 

032  =(-j2/ii!:a)yi(ao)  .  (58) 


=  f^VoW -  4  -  ^ 

V  2o  ^  o^  4 


X  2  \ 

J'i(W+[^Jr2(Ai)+  ^ kofc) .  (60) 


041  =  i-i2^ka)Ji{(xb)  , 
a42=i-i2nka)Yi{ab)  , 


«43  = 


yo03i>)+U^-^-^  Ji(eb)+^-^y2{pb)+  ^  h(fib)  .  (63) 


yo(P6)+Lt2_^_M 

r  4 


In  the  medium-wavenumber  region. 


CD  1,1  fl) 

—  <\k\<—  , 


and  the  matrix  coefficients  are 
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f  2  \  1  /X  \  f 

an=  (A  +  2/i)  —  -Xk-  /o(CM)+  —  /i(ow)  +  (A  +  2/i)  l2{aa)  .  (66) 

L  I  2  J  J  V  ;  I,  2  J 

f  2\  1 

fll2  =  (A+2/i)  -A/:^  ^q(om)  +  [ - \K\(cca)  +  {X  +  2f^)  — —  Al2{ow)  ,  (67) 

I  2  J  J  y  a  J  [2  ) 

fll3  =  {-inkp)jQiPa)  +  U^kP)J2ipa)  .  (68) 

fl,4  =  i-mkP)YQ{Pa)  +  mp)Y2iPa)  .  (69) 

/2'\  1  (^2^ 

fl2i=  (A+2/i)  ^  -Xk-  /o(a6)+  ^  /,(at)  +  (A  +  2/i)  ^  /2(afc) ,  (70) 

I  2  J  J  y  b  J  I  2  J 

022  =  iX  +  2n)  ^  -Xk^  Koiab)  +  \-^]Kl{ab)  +  iX  +  2^)  —  K2(ab)  ,  (71) 

I  2  j  J  \  b  J  {2  ) 

a2i=i-ifikP)Jo{Pb)  +  (i^ikP)J2iPb)  ,  (72) 

024  =  i-mkP)YQiPb)  +  mp)Y2{pb)  .  (73) 

^31  =(/2/iA:a)/i(CM)  ,  (74) 

^32  =(-/2;i/:a)A:i(afl)  ,  (75) 


041  =  ii2^ka)Iiiab)  - 
042  ={-i2nka)Kx{ab)  , 


J'i(A>)  +  (^]»’2(W+  ^  1'3(Ai)  .  (77) 


fl43=(ff)7o(/JW  +  U^-^-^ko3fc)  +  (-^j/2(W+  ^  73(W  ,  (80) 
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044 


=  (f>oW. 


F  4 


hW  •  (81) 


In  the  high-wavenumber  region, 
0)  0) 

and  the  matrix  coefficients  are 


flit  = 


fll2  = 


• 

(Af2/i) 

ce 

0 

-Ait2 

- 

1  2  J 

- 

- 

• 

a+2/i) 

a 

-Aik2 

<  2  , 

/o(afl)  +  |  —  j/i(aa)  +  U  +  2u)  —  hiouj)  , 

a  J  [2  J 


K 


KQ(aa)  + 


(=?) 


i:i(aa)-i-(A-i-2/i« 


v2y 


au=mP)Ko(pa)  +  mp)K2m  , 


fl2i  = 


^22  = 


(A +2^) 


V  2  y 


a 


2\ 


(A  +  2;/)  ^ 


-Ait^ 


-AA^ 


f  Xa\ 

I  b  ) 


Ioiab)+\^\lliab)  +  a  +  2nl\ 


a 


2^ 


hiab)  , 


\Ko(ab)+ 


(t) 


a 


2\ 


\Kiiab)  +  a+2ul^^ 


a23={-inkP)lQ{Pb)  +  {-inkp)l2iPb)  , 

a2A=mp)Koipb)+mp)K2m , 

fl3l  =  (i2nka)Ii{aa)  , 
a22=(-i2^lka)Klicca)  , 


^33  =f^ 


^  a_ti 
‘  2-"  4 


M- 


fl 


«34 


fl^  4 


/i03«)+(^]/2W' 


KiiPd)^ 


r^p2^ 


041  =  {i2^ka)Iiicxb)  » 


(82) 


(83) 


#:2(afl)  ,  (84) 

(85) 

(86) 


(87) 


\K2i0cb)  ,  (88) 

(89) 

(90) 

(91) 

(92) 


li(pa)  ,  (93) 


K^(Pa),  (94) 
(95) 
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a^2=(ri2iika)K\(ab)  , 


«43  =f^  koW  + 


F  4 


j'lW  +  (|f]/2(W' 


(96) 

Upb)  ,  (97) 


and 


044  =(^]A:o03fc)+^/it2  -^+^j«i(^W+(^)Ar2W 


*:3(/36)  .  (98) 


Finally,  the  displacements  of  the  forced  system  can  be  determined  by  use  of  the 
constants  Cj,  C2,  C3,  and  C4  in  the  low-wavenumber  region  as 


ikxJoM 


and 


Ur  =  [-C\CU\  (ar)  -  C2Ccyi  (ar)  -  C-^ikJi  {fir)  -  C4/7:yi  (y3r)]  e‘^e 


^2=1  C\ikJQ{ar)  +  C2ikYQ{ar) 


(f  ]-/o  m  +  0r)  -  (^|]72(/3'-) 


JkxJast  . 
”  e  ; 


+ 


+  C4  [|]l'o(P'-)  +  (^ijj’l(Pr)-(^|y203'-)  } 

in  the  medium-wavenumber  region  as 

Ur  =  [Cia/i(ar)  -  C2aKi{ar)  -  C^iUiipr)  -  C4ikYiipr)] 
and 

«2  ={  CiiklQiar)  +  C2ikKQiar) 


+  C3 


+  C4 

and  in  the  high-wavenumber  region  as 


(f  j-'oOS'-) + (7)^1  (M  -  (f  )-'2(w] 
(^1  jKo(Pr) jKi(Pr)  -  [|y2(P'') 


JkxJwt 
e  e 


(99) 


(100) 


(101) 


(102) 
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Ur  =  [CiOdiiar)  -  C2aKiiar)  -  C^ikliipr) -  C4ikKiiPr)] 
and 

“z  ={  CiiklQiar)  +  C2ikKQiar) 


+  ^3 


+  C4 


(f  ]/o(P'-)  +  (;)/l  (Pr) + [|]/20'-)] 

{^^00'-) + -  (f  ]jf2(P'-)] 


Equations  (99),  (101),  and  (103)  can  be  written  as 


\  ^ikxjox 
]  e  e 


Ur  =  Urir)e^e‘^  , 

and  equations  (100),  (102),  and  (104)  can  be  expressed  as 
=  V^{r)e^e'^  . 


(103) 


(104) 


(105) 


(106) 


3.  MODEL  VALIDATION 


The  new  thick  shell  model  is  validated  by  comparing  it  to  a  thin  shell  problem.  The 
thickness  of  the  thick  shell  model  is  set  to  a  very  small  value  (relative  to  the  shell  radius)  to 
ensure  such  an  accurate  comparison.  The  governing  partial  differential  equations  for  an 
axisymmetric,  isotropic  thin  shell  are  (Hull,  1994) 


ph 


d^Urizj) 

dt^ 


h^Ej^  d^Uriz,t) 

12  ~dF~ 


-^Ur{z,t) 


a  az 


(108) 
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where  is  the  effective  modulus  of  elasticity,  u  is  Poisson's  ratio,  a  is  the  mean  radius  of  the 
shell,  and  h  is  its  thickness.  The  relationship  between  the  effective  modulus  and  Young's 
modulus  is  given  by 


where  E  is  Young's  modulus.  The  definition  of  the  Lamd  constants  (from  section  2)  in  terms 
of  material  properties  (Young's  modulus  £  and  Poisson's  ratio  u)  are  (Mirsky  and  Herrmann, 
1958) 


A  = 


vE 

(l+o)(l-2t;) 


(110) 


and 


£ 

2(1 +  u)  ■ 


(111) 


Equations  (107)  and  (108)  are  solved  as  an  infinite  length  system  by  setting  the 
displacements  and  forcing  functions  at  a  specific  wavenumber  and  frequency,  resulting  in  the 
following  equations: 


Ur(x,t)  =  Ure‘^e‘^  ,  (112) 

=  ,  (113) 

f{x,t)  =  Fe^e'^  ,  (114) 

and 

p(x,t)  =  pA'®' .  (115) 


Inserting  the  values  and  derivatives  of  equations  (112)-(115)  into  equations  (107)  and  (108) 
produces  a  two-by-two  system  of  linear  algebraic  equations: 


’^1 

p2\  hlA^r]  l^J 

where 


(116) 
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1  =  hEjlP"  -  hpoP"  , 

(117) 

,  -ikhvE, 

bi2  = 

a 

(118) 

,  -ikhvEr 

hi= - ^  • 

(119) 

a 

.  ^  2  hE,  h^E,k^ 

b22=ph(0  — . 

(120) 

Note  that  both  and  V,  are  constant  with  respect  to  r  for  the  thin  shell  equations  of  motion. 

Figures  3-6  show  the  transfer  functions  of  displacements  divided  by  forcing  functions 
for  a  shell  with  a  mean  radius  of  0.1  m,  thickness  of  0.001  m,  Poisson’s  ratio  of  0.4 
(dimensionless).  Young's  modulus  of  1  x  10^  N/m^,  density  of  2000  kg/m^,  and  frequency  of 
50  Hz.  The  ratio  of  thickness  to  radius  {hid)  is  1/100.  The  corresponding  dilatational  wave 
speed  is  1035.1  m/s  and  the  corresponding  shear  wave  speed  is  422.6  m/s.  In  these  figures, 
the  solid  line  represents  the  thick  shell  model  and  the  X's  denote  the  thin  shell  model.  Figure  3 
plots  the  magnitude  of  the  axial  displacement  divided  by  applied  axial  stress  versus 
wavenumber,  and  figure  4  plots  the  magnitude  of  the  radial  displacement  divided  by  applied 
axial  stress  versus  wavenumber.  In  both  figures,  the  applied  radial  pressure  on  the  shell  was 
set  equal  to  zero.  Figure  5  is  a  plot  of  the  magnitude  of  the  axial  displacement  divided  by 
applied  radial  pressure  versus  wavenumber,  and  figure  6  is  a  plot  of  the  magnitude  of  the  radial 
displacement  divided  by  applied  radial  pressure  versus  wavenumber.  In  both  figures,  the 
applied  axial  stress  on  the  shell  was  set  equal  to  zero.  All  the  figures  show  the  thick  shell 
model  to  be  in  agreement  with  the  thin  shell  model.  Because  the  figures  ate  symmetric  in 
wavenumber,  the  negative  wavenumbers  are  not  shown.  The  thick  shell  forcing  functions 
were  applied  to  the  exterior  of  the  thick  shell  model.  Identical  results  would  be  obtained  by 
application  of  the  loads  to  the  interior  of  the  thick  shell  because  the  shell  thickness  in  this 
example  is  set  to  such  a  small  value. 
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0  0.5  1.0  1.5  2.0 

Wavenumber,  k  (rad/m) 

Figure  3.  Magnitude  of  Axial  Displacement  Divided  by  Applied  Axial  Stress: 
Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model, /=  50  Hz,  h/a  =  1/100 


Figure  4.  Magnitude  of  Radial  Displacement  Divided  by  Applied  Axial  Stress: 
Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/=  50  Hz,  h/a  =  1/100 
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Wavenumber,  k  (rad/m) 

Figure  5.  Magnitude  of  Axial  Displacement  Divided  by  Applied  Radial  Pressure: 
Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/  =  50  Hz,  h/a=  II 100 
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Thick  Shell  Model 


Thin  Shell  Model 
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Figure  6.  Magnitude  of  Radial  Displacement  Divided  by  Applied  Radial  Pressure: 
Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model,/=  50  Hz,  hla  =  1/100 
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4.  ANALYSIS 

The  thick  shell  model  is  now  investigated  for  a  shell  whose  thickness  to  mean  radius 
ratio  [(2/i/(a+/j)]  is  1.  For  this  example,  the  mean  radius  [(fl+Z?)/2]  is  0.1  m,  the  inner  radius 
(a)  is  0.05  m,  the  outer  radius  {b)  is  0.15  m,  and  the  shell  thickness  (/i)  is  0.1  m.  Other 
properties  of  the  shell  are  a  Poisson's  ratio  of  0.4  (dimensionless),  a  Young's  modulus  of  1  x 
10^  N/m^,  and  a  density  of  2000  kg/m^.  The  forcing  functions  were  applied  to  the  exterior  of 
the  shell  at  two  frequencies:  50  and  500  Hz. 

In  figures  7-14,  the  thick  shell  model  is  compared  to  the  thin  shell  model.  The  thick 
shell  model  is  shown  at  the  mean  radius  value  of  r  =  {a+b)l2  and  denoted  using  a  solid  line, 
while  the  thin  shell  model  is  shown  with  a  dashed  line.  Figure  7  is  a  comparison  of  the 
magnitude  of  the  axial  displacement  divided  by  applied  axial  sU'ess  versus  wavenumber,  and 
figure  8  is  a  comparison  of  the  magnitude  of  the  radial  displacement  divided  by  applied  axial 
stress  versus  wavenumber.  In  both  figures,  the  applied  radial  pressure  on  the  shell  was  set 
equal  to  zero.  Figure  9  plots  the  rr.ignitude  of  the  axial  displacement  divided  by  applied  radial 
pressure  versus  wavenumber,  and  figure  10  plots  the  magnitude  of  the  radial  displacement 
divided  by  applied  radial  pressure  versus  wavenumber.  In  both  figures,  the  applied  axial  stress 
on  the  shell  was  set  equal  to  zero.  The  analysis  in  figures  7-10  was  computed  at  50  Hz.  In 
most  of  the  wavenumbers  shown,  the  thin  shell  model  underpredicts  the  system  response  by  2 
to  8  dB.  Figure  1 1  compares  the  magnitude  of  the  axial  displacement  divided  by  applied  axial 
stress  versus  wavenumber,  and  figure  12  compares  the  magnitude  of  the  radial  displacement 
divided  by  applied  axial  stress  versus  wavenumber.  In  both  figures,  the  applied  radial  pressure 
on  the  shell  was  set  equal  to  zero.  Figure  13  is  a  plot  of  the  magnitude  of  the  axial 
displacement  divided  by  applied  radial  pressure  versus  wavenumber,  and  figure  14  is  a  plot  of 
the  magnitude  of  the  radial  displacement  divided  by  applied  radial  pressure  versus 
wavenumber.  In  both  figures,  the  applied  axial  stress  on  the  shell  was  set  equal  to  zero.  The 
analysis  in  figures  11-14  was  computed  at  500  Hz.  Figures  11-14  show  behavior  at  low 
wavenumbers  similar  to  that  shown  in  figures  7-10  where  the  thin  shell  theory  usually 
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underpredicts  the  response.  There  is  a  crossover  of  the  thick  shell  and  the  thin  shell  models  at 
higher  wavenumbers  in  figures  1 1  and  12. 

In  figures  15-22,  the  thick  shell  model  response  at  the  exterior  of  the  shell  (r  =  b)  is 
compared  to  the  thick  shell  model  response  at  the  interior  of  the  shell  (r  =  a).  The  solid  line 
denotes  the  response  at  r  s  a  and  the  dashed  line  at  r  =  b.  Figure  15  compares  the  magnitude  of 
the  axial  displacement  divided  by  applied  axial  su-ess  versus  wavenumber,  and  figure  16 
compares  the  magnitude  of  the  radial  displacement  divided  by  applied  axial  stress  versus 
wavenumber.  Figure  17  is  a  plot  of  the  magnitude  of  the  axial  displacement  divided  by  applied 
radial  pressure  versus  wavenumber,  and  Hgure  18  is  a  plot  of  the  magnitude  of  the  radial 
displacement  divided  by  applied  radial  pressure  versus  wavenumber.  The  analysis  in  figures 
15-18  was  computed  at  50  Hz.  These  figures  correspond  to  figures  7-10,  except  the  values 
now  shown  are  the  transfer  functions  at  the  interior  and  the  exterior  of  the  shell.  For  thin  shell 
theory,  these  two  transfer  functions  are  identical  (and  also  equal  to  the  value  at  the  mean  radius 
of  the  shell).  Note  that  in  figure  16  these  two  transfer  functions  are  significantly  different. 
Figure  19  is  a  comparison  of  the  magnitude  of  the  axial  displacement  divided  by  applied  axial 
stress  versus  wavenumber,  and  figure  20  is  a  comparison  of  the  magnitude  of  the  radial 
displacement  divided  by  applied  axial  stress  versus  wavenumber.  Figure  21  plots  the 
magnitude  of  the  axial  displacement  divided  by  applied  radial  pressure  versus  wavenumber, 
and  figure  22  plots  the  magnitude  of  the  radial  displacement  divided  by  applied  radial  pressure 
versus  wavenumber.  The  analysis  in  figures  19-22  was  computed  at  500  Hz.  These  figures 
correspond  to  figures  11-14,  except  the  values  now  shown  are  the  transfer  functions  at  the 
interior  and  the  exterior  of  the  shell.  Note  that,  except  for  figure  22,  the  transfer  function  at  the 
interior  of  the  shell  compared  with  that  at  the  exterior  of  the  shell  is  significantly  different 
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Figure  7.  Magnitude  of  Axial  Displacement  Divided  by  Applied  Axial  Stress: 
Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model, /=  50  Hz,  2hJia+b)  =  1 


Figure  8.  Magnitude  of  Radial  Displacement  Divided  by  Applied  Axial  Stress: 
Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model, /=  50  Hz,  2hl(,a+b)  =  1 
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Figure  9.  Magnitude  of  Axial  Displacement  Divided  by  Applied  Radial  Pressure: 
Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model, /=  50  Hz,  2h/{a+b)  =  1 


Figure  10.  Magnitude  of  Radial  Displacement  Divided  by  Applied  Radial  Pressure: 
Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model, /=  50  Hz,  7hJia+b)  =  1 
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Figure  1 1.  Magnitude  of  Axial  Displacement  Divided  by  Applied  Axial  Stress: 
Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model, /=  500  Hz,  2hJ(a-¥b)  =  1 


Figure  12.  Magnitude  of  Radial  Displacement  Divided  by  Applied  Axial  Stress: 
Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model, /=  500  Hz,  2hlia+b)  =  1 
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Figure  13.  Magnitude  of  Axial  Displacement  Divided  by  Applied  Radial  Pressure: 
Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model, /=  500  Hz,  2hl{a+b)  =  1 


Figure  14.  Magnitude  of  Radial  Displacement  Divided  by  Applied  Radial  Pressure: 
Comparison  of  Thick  Shell  Model  to  Thin  Shell  Model, /=  500  Hz,  2h/(a+fc)  =  1 
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Figure  15.  Magnitude  of  Axial  Displacement  Divided  by  Applied  Axial  Stress: 
Comparison  of  Transfer  Function  at  r  =  a  to  r  =  fc,/=  50  Hz,  2hlia+b)  =  1 


Figure  16.  Magnitude  of  Radial  Displacement  Divided  by  Applied  Axial  Stress: 
Comparison  of  Transfer  Function  at  r  =  a  to  r  =  b,f=  50  Hz,  2hHja-tb)  =  1 
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Figure  17.  Magnitude  of  Axial  Displacement  Divided  by  Applied  Radial  Pressure: 
Comparison  of  Transfer  Function  at  r  =  a  to  r  =  b,f=  50  Hz,  2hJ(a+b)  =  1 


Figure  18.  Magnitude  of  Radial  Displacement  Divided  by  Applied  Radial  Pressure: 
Comparison  of  Transfer  Function  at  r  =  a  to  r  =  fc,/=  50  Hz,  2h/(a+b)  =  1 
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Figure  19.  Magnitude  of  Axial  Displacement  Divided  by  Applied  Axial  Stress: 
Comparison  of  Transfer  Function  at  r  =  a  to  r  =  i),/=  500  Hz,  2h/ia+b)  =  1 


Figure  20.  Magnitude  of  Radial  Displacement  Divided  by  Applied  Axial  Stress: 
Comparison  of  Transfer  Function  at  r  =  a  to  r  =  b,f=  500  Hz,  2hlia+b)  =  1 
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Figure  21.  Magnitude  of  Axial  Displacement  Divided  by  Applied  Radial  Pressure: 
Comparison  of  Transfer  Function  at  r  =  a  to  r  =  500  Hz,  2h/ia+b)  =  1 


Figure  22.  Magnitude  of  Radial  Displacement  Divided  by  Applied  Radial  Pressure: 
Comparison  of  Transfer  Function  at  r  =  a  to  r  =  b,f=  500  Hz,  2hlia+b)  =  1 
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5.  CONCLUSIONS 

A  model  of  an  infinite  length,  axisymmetric,  isotropic,  forced  thick  shell  response 
based  on  Bessel  function  solutions  has  been  presented.  This  new  model  allows  for 
displacement  variations  across  the  shell  wall.  The  case  of  the  mean  radius  of  the  shell  equal  to 
the  shell  thickness  has  been  analyzed  and  compared  to  thin  shell  equations  of  motion.  It  is 
shown  that,  in  general,  the  thin  shell  equations  of  motion  tend  to  underpredict  thick  shell 
response. 
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APPENDIX  -  DERIVATIVES  OF  THE  BESSEL  FUNCTIONS 

The  following  equations  are  derivatives  of  the  Bessel  functions  used  in  section  2: 


^[yo(Ofr)]  =  -a/j(ar)  , 

(A-1) 

~[yi(ar)]  =  ^[jQiar) -  J2{ar)]  . 

(A-2) 

^[y2(Of»')]  =  -j[^l(ar)-y3(ar)]  . 

(A-3) 

■^[YQ{ar)]  =  -aYx{ar)  , 
dr 

(A-4) 

j^[Yx{ar)]  =  SL[YQ{ar)-Y2{ar)]. 

(A-5) 

jl[Y2{ar)]  =  ^[Y^{ar)-Y^{ar)]  , 

(A-6) 

■^[KQ{ar)]  =  -aKi{ar)  , 

(A-7) 

±[Kx{ar)]^^[K2{ar)-\-KQ{ar)]  . 

(A-8) 

^[/i:2(ar)]  =  :^[/i:3(ar)  +  /(:i(ar)]  , 

(A-9) 

^[lQ{ar)]  =  -aIx{ar)  , 

(A-10) 

^[/l(ar)]  =  -j[/2(ar)  +  /o(ar)]  . 

(A-11) 

^[/2(ar)]  =  •|[/3(ar)  +  /i(ar)]  . 

(A-12) 

A-l/A-2 
Reverse  Blank 
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